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Landau’s phenomenological theory of Fermi liquids is a fundamental paradigm in many-body 
physics that has been remarkably successful in explaining the properties of a wide range of interacting 
fermion systems, such as liquid helium-3, nuclear matter, and electrons in metals. The d-dimensional 
boundaries of (d -|- l)-dimensional topological phases of matter such as quantum Hall systems and 
topological insulators provide new types of many-fermion systems that are topologically distinct 
from conventional d-dimensional many-fermion systems. We construct a phenomenological Landau 
theory for the two-dimensional helical Fermi liquid found on the surface of a three-dimensional time- 
reversal invariant topological insulator. In the presence of rotation symmetry, interactions between 
quasiparticles are described by ten independent Landau parameters per angular momentum chan¬ 
nel, by contrast with the two (symmetric and antisymmetric) Landau parameters for a conventional 
spin-degenerate Fermi liquid. We then project quasiparticle states onto the Fermi surface and ob¬ 
tain an effectively spinless, projected Landau theory with a single projected Landau parameter per 
angular momentum channel that captures the spin-momentum locking or nontrivial Berry phase 
of the Fermi surface. As a result of this nontrivial Berry phase, projection to the Fermi surface 
can increase or lower the angular momentum of the quasiparticle interactions. We derive equilib¬ 
rium properties, criteria for Fermi surface instabilities, and collective mode dispersions in terms of 
the projected Landau parameters. We briefly discuss experimental means of measuring projected 
Landau parameters. 

PACS numbers: 71.10.-w, 71.10.Ay, 71.70.Ej, 73.20.-r 


The Landau theory of Fermi liquids (FL) [ 1 ], or FL 
theory for short, is the cornerstone of our understand¬ 
ing of weakly correlated, gapless Fermi systems at low 
temperatures, such as ^He atoms in the normal liquid 
state and itinerant electrons in metals. FL theory ex¬ 
plains the puzzling observation that despite strong inter¬ 
actions between the constituent fermions, many Fermi 
systems behave essentially as free Fermi gases, except for 
the renormalization of their physical properties which is 
captured by dimensionless quantities known as Landau 
parameters. These Landau parameters describe how the 
elementary excitations of the FL—the quasiparticles and 
quasiholes—interact with one another. 

Topological insulators [ 2 ] provide new types of gapless 
Fermi systems: topological surface/edge states. In the 
absence of interparticle interactions, electrons propagat¬ 
ing on the edge of a two-dimensional ( 2 D) topological 
insulator [ 3 ] form a ID helical Fermi gas [ 4 ]. In the pres¬ 
ence of interactions, the ID helical Fermi gas becomes 
a ID helical Luttinger liquid [ 5 ] with no sharply defined 
Fermi points. In 3 D topological insulators, surface elec¬ 
trons form a 2 D helical Fermi gas [6] , which is expected to 
evolve adiabatically into a 2 D helical FL in the presence 
of electron-electron interactions. 

This paper presents a FL theory for the interacting 
2 D surface states of the 3 D topological insulator. To 
our knowledge, such a helical FL theory has been miss¬ 
ing in the literature despite the recent surge of interest 


in the effects of electron-electron interactions in topo¬ 
logical insulators [ 7 ]. In the spirit of standard FL the¬ 
ory [1], we focus on systems with a discrete time-reversal 
symmetry, the protecting symmetry of topological in¬ 
sulators, as well as continuous translation and spatial 
rotation symmetries. We further consider the simplest 
case of a single surface Fermi surface—denoted simply 
as the Fermi surface in the following—which by rota¬ 
tion symmetry must be circular. This does not apply 
to certain topological insulators whose Fermi surface is 
strongly anisotropic, such as Bi2Te3 with 0 . 67 % Sn dop¬ 
ing [8] where there are large hexagonal warping effects 
due to the rhombohedral crystal structure of the bulk ma¬ 
terial [ 9 ]. However, in several other topological insulators 
such as Bi2Se3 [10], Bi2Te2Se [11], SbxBi2_a;Se2Te [11], 
Bi1.5Sbo.5Te1.7Se1,3 [ 12 ], Tli_j,Bii+3;Se2-5 [ 13 ], strained 
a-Sn on InSb(OOl) [ 14 ], and strained HgTe [ 15 ], the Fermi 
surface as observed in angle-resolved photoemission spec¬ 
troscopy (ARPES) is very nearly circular. However, 
due to spin-momentum locking in the topological sur¬ 
face states [6]— a consequence of strong spin-orbit cou¬ 
pling, rotation symmetry in a helical FL must neces¬ 
sarily involve spin degrees of freedom, which leads to 
a theory rather different from that of the conventional 
spin-degenerate FL. Moreover, the existence of a single 
nondegenerate Fermi surface—a consequence of the topo¬ 
logical character of the bulk—eventually leads, via the 
application of the general principles of FL theory, to an 
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effectively spinless FL theory. The physical properties 
of the resulting helical FL are nevertheless distinct from 
those of a truly spinless FL, due to a nontrivial mapping 
between physical, spinful quasiparticles, and the effec¬ 
tive, spinless quasiparticles. For the same reason, our 
helical FL theory is also qualitatively different from re¬ 
cently constructed FL theories of non-topological spin- 
orbit coupled systems such as the Rashba 2D electron 
gas [16] and 3D spin-orbit coupled metals [17], which are 
characterized by two (spin-split) Fermi surfaces. 

FL theory views the many-fermion system as a gas of 
elementary excitations above the ground state, the quasi¬ 
particles. Because translation symmetry is assumed, the 
momentum p = {px,Py) of the quasiparticles is well- 
defined and a configuration of quasiparticles is specified 
by a distribution function Up. In a conventional FL, spin 
is conserved and the distribution function is diagonal in 
spin space Upa = {clp^Cpffj, where c],g. {Cpff} is a creation 
(annihilation) operator for a fermion with momentum p 
and spin tr =t,4-, but in systems with spin-orbit coupling 
such as the helical FL the distribution function is gen¬ 
erally a matrix in spin space, = (cj^^Cp^) [17] . The 
central quantity in FL theory is the energy 5E of the 
gas of interacting quasiparticles relative to the ground- 
state energy, expressed as a functional of the deviation 
Sup^ = Up^ — np'^°‘^ of the distribution function from its 
value in the ground state, 

SE[ 6 np]= J dpha 0 {p)Snp^ 

+ lj dpAp' Va0-,js{p,p')Sn'i^^Snl^, ( 1 ) 

where (working in units such that h = 1 ) 

h{p) = vpz ■ (cr X p), (2) 

is the single-particle Dirac Hamiltonian of the topolog¬ 
ical surface state [2] with vp the Fermi velocity [18], 
yai 3 -js{p,p') is a reduced two-body interaction that de¬ 
pends only on the unit vector p = p/\p\ parameterizing 
the Fermi surface, and we denote the integration mea¬ 
sure by / dp = J . The form of Eq. (1) can be 
obtained from a generic, translationally invariant inter¬ 
action Vai 3 -js{k, q) by requiring that all fermionic mo¬ 
menta lie on the Fermi surface [19]. 

Our first goal is to derive the most general form of 
the two-body interaction Vh/ 3 ;-y 5 (p, p') consistent with the 
general principles of quantum mechanics and the symme¬ 
tries of the problem. This goal is most easily achieved by 


expanding the two-body interaction as 

3 oo 

Vap--yS{p,p') = Y ( 3 ) 

fj,p=o oo 

where p = (cos0p, sin0p), 1 , 1 ' are angular momentum 
quantum numbers, and the set of four 2x2 Hermitian 
matrices = (1 , <t) where 1 denotes the identity matrix 
allows us to construct the quasiparticle charge Spp and 
spin 6 Sp densities {i = x, y, z), 

5pp = a%5nf = 5^p5nf, 6 s], = \alp5nf. (4) 

Upon substituting Eq. (3) in Eq. (1), one obtains three 
classes of terms: charge-charge interactions proportional 
to Vqq , spin-spin interactions proportional to Vfj , and 
spin-charge interactions proportional to VqI = V^qK 
Time-reversal symmetry implies that the angular mo¬ 
menta I and V must differ by an even integer for charge- 
charge and spin-spin interactions and by an odd integer 
for spin-charge interactions [19] . 

The main difference between a conventional FL and a 
spin-orbit coupled FL such as the helical FL lies in the 
consequences of rotation symmetry. The single-particle 
Hamiltonian (2) is neither invariant under a spatial ro¬ 
tation nor under a spin rotation, but is invariant under 
a simultaneous rotation of spatial and spin coordinates: 
[Jz,h{p)] = 0, where Jz = —*gf- + is the total (or¬ 
bital plus spin) angular momentum in the z direction. 
Requiring that the interaction term in Eq. (1) be also 
invariant under such rotations, we find that it can be 
written as the sum of three terms 6 Vcc, SVac, and 514,,, 
where [19] 

OO 

1=0 


dp dp' fff cos I0pp> SppSpp ,, (5) 


is the charge-charge interaction. 


514c = 


OO 


Z =0 * 


dp dp' 


X [ifY cos Wpp, -I- sin Wpp,)Sppp' ■ Ssp> 

+ cosWpp’ + smWppt)6ppp' X Sspt 


( 6 ) 


is the spin-charge interaction, and 
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= 5 X! / 3p<l/{cosZ0pp/ (/r’^(<5Sp(5Sp, + Ss^Ss'^,) + fr'‘^5sp8Sp^ + sinZ6lpp/i5sp x dsp, 

1=0 

+ cosldpp>(^fi''’^ [{p ■ Ssp) {p' X 6sp,) + (p X (5sp) (p' • (5sp/)] + HP ' iP'' '^V) “ (P x Ssp) (p' x (5v)]) }i 

( 7 ) 


is the spin-spin interaction. We denote by 0pp' = Op/ — 6 p 
the relative angle between p and p', and write a x b = 
z ■ {a X b) for the cross product of two in-plane vectors. 

Equations (5)-(7), the first main result of this work, 
represent the most general short-range two-body inter¬ 
action in a helical FL consistent with translation, ro¬ 
tation, and time-reversal symmetries. The interaction 
is specified by ten real Landau parameters for each 
value of the relative angular momentum I = 0 , 1 , 2 ,...: 
one charge-charge parameter /“, four spin-charge pa¬ 
rameters . • ■,and five spin-spin parameters 

..., This stands in contrast to the two Lan¬ 

dau parameters ff (spin symmetric) and /“ (spin anti¬ 
symmetric) in a conventional FL [1], which would cor¬ 
respond to // = fr and ft = yr" = 
the absence of spin-orbit coupling. In particular, spin- 
orbit coupling allows for a nonzero spin-charge inter¬ 
action ( 6 ) which would be forbidden by separate spa¬ 
tial and spin rotation symmetries in a conventional FL. 
The spin-spin interaction (7) also exhibits novel features: 
fi^'^ 7 ^ general, which corresponds to an XXZ 

interaction with Ising anisotropy rather than the conven¬ 
tional SU (2)-symmetric Heisenberg interaction; is 

a Dzyaloshinskii-Moriya interaction; and are 

anisotropic spin-spin interactions similar to those found 
in compass models [ 20 ], but with a continuous rather 
than discrete spin-orbit rotation symmetry. 

While Eq. (5)-(7) in conjunction with Eq. (1) cor¬ 
rectly describe the helical FL, in the spirit of FL the¬ 
ory one can go one step further and only retain elec¬ 
tron states on the Fermi surface. Because of the strong 
spin-orbit coupling present in the Dirac Hamiltonian (2), 
such electrons are annihilated by the operator ijjp± = 
-^(ie“*®pcp-|- ± Cp^), where positive (-I-) helicity corre¬ 
sponds to a positive Fermi energy ep > 0 above the Dirac 
point, and negative (—) helicity corresponds to a nega¬ 
tive Fermi energy ep < 0. Inverting this relation, one can 
express the spin eigenoperators Cp^ in terms of the helic¬ 
ity eigenoperators 'tj;p± as Cp^ = y^^iippp + tfp-) and 
Cpi = ■^('!/ip_|_ — ifp-). Choosing ep > 0 for definiteness, 
the Fermi surface consists exclusively of electron states 
of positive helicity, such that one may wish to drop the 
negative helicity eigenoperators ifp- entirely from these 
expressions for Cp-j- and Cpj,. Applying this procedure to 
Eq. ( 1 ) yields a Landau functional for an effectively spin¬ 


less FL theory. 



+ dpdp'fI COS lOpp^SupSupr, ( 8 ) 

1^0 

where e® = vp\p\ is the dispersion relation of positive he¬ 
licity quasiparticles, Sfip = fip—np with fip = (';/'],+V'p+) 
is the distribution function for these quasiparticles, and 
fi are effectively spinless, projected Landau parameters 
related to the ten unprojected Landau parameters previ¬ 
ously discussed by 

n _ rCC rSC,3 1 pSS,5 

Ji — Ji ~ h ~ ih 

+ W!y-f!y + f!y + f!iy ( 9 ) 

for I = 0 , 1 , 2 ,..., with the definition ^ = 0 . 

The quasiparticle charge and spin densities (4) are given 
in terms of Sfip by 

6pp = Snp, 5sp = \eijPj5fip, i = x,y, Ssp = 0, (10) 

where the last two equalities express spin-momentum 
locking in the xy plane. Equations (8)-(10), together 
with the definitions of the unprojected Landau param¬ 
eters in Eq. (5)-(7), are the second main result of this 
work. 

Before deriving the physical properties of the helical 
FL from the projected Landau functional ( 8 ), we pause 
to discuss a number of interesting features of the re¬ 
lationship (9) between projected and unprojected Lan¬ 
dau parameters. The unprojected Landau parameters 
and do not enter the projected interac¬ 
tion because spin and momentum are perpendicular on 
the Fermi surface (p • Ssp = 0) due to spin-momentum 
locking. The parameter does not enter either be¬ 
cause it produces a projected interaction that is odd un¬ 
der pop', which is inconsistent with particle indis- 
tinguishability. The last term on the right-hand side of 
Eq. (9) shows that projection to the Fermi surface can 
effectively raise or lower the angular momentum of the 
unprojected interaction. For example, for I = 1 one has 

_ rCC rSC,3 1 rSS,5 

— Jl Jl 4/1 

+ l(/o“’'-/o*’'+/r'+/r'), ( 11 ) 
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that is, an isotropic, s-wave {I = 0 ) microscopic interac¬ 
tion can produce an anisotropic, p-wave (/ = 1 ) effective 
interaction in the projected theory. This can be seen as 
the particle-hole counterpart to the effective p-wave in¬ 
teraction in the Bardeen-Cooper-Schrieffer (BCS) chan¬ 
nel produced on the doped surface of a 3D topological 
insulator by a microscopic s-wave BCS interaction [21]. 

As in standard FL theory, many physical properties of 
the helical FL can be derived from the projected Landau 
functional ( 8 ). The simplest property is Luttinger’s theo¬ 
rem [22], i.e., the relation pi?- = y/inn between Fermi mo¬ 
mentum pp and total density n of quasiparticles, which 
is also equal to the total density of electrons (defining a 
system with pp = 0 as the vacuum). That Luttinger’s 
theorem holds in its original form despite the presence of 
strong spin-orbit coupling is a consequence of the exis¬ 
tence of a single helical Fermi surface, which is only possi¬ 
ble on the surface of a 3D topological phase. Interactions 
in topologically trivial spin-orbit coupled systems such 
as the Rashba 2D electron gas can individually renor¬ 
malize the Fermi momenta of the two spin-split Fermi 
surfaces [16]. Other equililibrium properties of the heli¬ 
cal FL can be calculated from the quasiparticle energy 
Cp, defined as the functional derivative of the Landau 
functional with respect to the distribution function. 


SE 

5fii 


= 


E 

1=0 • 


dp'/; cos WppiSupi. ( 12 ) 


From Eq. (12) one can follow the standard FL ap¬ 
proach [19] to derive the electronic specific heat coeffi¬ 
cient 7 = Cy/T and electronic compressibility k of the 
helical FL at zero temperature. 


7 = \T^^k%p{ep), 


^ p(£f) 1 _ 

l + Fo 


(13) 


where we define dimensionless Landau parameters Fb = 
p{ep)fo and Fi = \p{ep)fi, Z = 1,2,3,..., with p{tp) = 
€pI2'i:v^ the density of states of the helical FL at the 
Fermi energy ep = vppp. The compressibility becomes 
negative for Fq < — 1 , signaling an instability towards 
phase separation [23]. Unlike in a standard FL, here this 
condition can be reached not only for attractive density- 
density interactions, but also as a result of spin-charge or 
even purely spin-spin interactions, given the relation (9) 
between the projected and unprojected Landau parame¬ 
ters. 

The renormalized Fermi velocity vp differs in gen¬ 
eral from the Fermi velocity of noninteracting electrons 
v%. This is similar in spirit to the renormalization of 
the quasiparticle mass in a standard FL. The deriva¬ 
tion of the latter relies on Galilean invariance, while in 
the helical FL, Galilean invariance is broken by spin- 
orbit coupling. However, adiabatic continuity still im¬ 
plies that the total flux of quasiparticles is equal to the 
total flux of electrons [1]. The latter is calculated from 


the quantum-mechanical velocity operator for electrons 
Ve = v^p{^ X cr) which, for momentum-independent mi¬ 
croscopic interactions [24], is the same as in the absence 
of interactions [25] : it is a function of the noninteracting 
Fermi velocity, rather than the renormalized one. The 
total quasiparticle flux is a function of the quasiparticle 
velocity Vqp = VpCp. Equating the two fluxes yields a 
relation between the two Fermi velocities [19], 

^ = 1 + Fi, (14) 

Vp 

which is the helical FL analog of the relation — = 
1 -|- iFf between renormalized m* and noninteracting 
m quasiparticle masses in a standard FL [1]. 

The spin susceptibility introduces some added sub¬ 
tleties: unlike in a standard FL, it is not, strictly speak¬ 
ing, a Fermi surface property. Indeed, it depends explic¬ 
itly on a high-energy cutoff A already in the noninter¬ 
acting limit [26, 27[. In a standard FL, one can always 
choose the spin quantization axis to be parallel to the 
applied magnetic field B, such that the quasiparticle en¬ 
ergy shift Sepa = ^gp^sBa due to Zeeman coupling {g is 
the g-factor, /rs is the Bohr magneton) is diagonal in the 
spin basis cr = ±1. The resulting change in occupation 
numbers is localized to the Fermi surface in the zero-field 
limit, causing the spin susceptibility to be a Fermi sur¬ 
face property. In the helical FL, there is no freedom to 
choose the spin quantization axis due to spin-momentum 
locking, and the Zeeman coupling contains off-diagonal 
terms in the helicity basis. The projected FL theory (8), 
which projects out negative helicity states, cannot take 
these off-diagonal terms into account and thus should not 
be expected to yield exact results for the spin suscepti¬ 
bility. Nevertheless, one can calculate the Fermi surface 
contribution to the spin susceptibility using (8) and com¬ 
pare it in the noninteracting limit to an exact calculation 
that takes both helicities into account. The spin suscep¬ 
tibility tensor Xij is found to be diagonal, with in-plane 
Xxx = Xyy out-of-plane Xzz components given by 

Xxx = y'^pjipiep)—^, Xzz = 0, (15) 

1 H" -Tl 

in the projected FL theory, and 
Xxx = y^plp{A), Xzz = jg^pl[p{A) - p{ep)], (16) 

for the noninteracting Dirac surface state, including both 
helicities [19[. Thus in the noninteracting limit, Eq. (15) 
and (16) agree in the formal limit of large Fermi en¬ 
ergy ep ^ A. By contrast with the spin susceptibility 
of the standard FL which is renormalized by the spin- 
antisymmetric I = 0 Landau parameter Fq, here it is 
renormalized by a / = 1 Landau parameter due to spin- 
momentum locking on the Fermi surface. 

Pomeranchuk instabilities [28[ are instabilities of the 
Fermi surface towards spontaneous, static distortions of 
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its shape. To study such instabilities in the helical FL, 
one characterizes distortions of the Fermi surface by an 
angle-dependent Fermi momentum, expanded in angular 
momentum components, 

OO 

Pf{0)-Pf= Y. 

1 — — 00 

where A-i = A* because pf{0) is real. Substituting this 
expression into the Landau functional ( 8 ), one finds that 
the energy 5E remains positive, and thus the helical FL 
stable, if and only if [19] 

Fi > - 1 , (18) 

for all I = 0,1, 2,... This is the same as Pomeranchuk’s 
original criterion in 2D, but applied this time to the pro¬ 
jected Landau parameters, which are nontrivial functions 
of the unprojected ones. It contains as special cases the 
instability towards phase separation, already seen, as well 
as an instability towards in-plane magnetic order [29] for 
Fi —>■ — 1, that is signaled by divergences of the in-plane 
spin susceptibility (15) and the renormalized Fermi veloc¬ 
ity (14). The latter divergence also accompanies the I = 1 
spin-symmetric Pomeranchuk instability of the standard 
FL [30] . The I = 2 instability is towards quadrupolar dis¬ 
tortions of the helical Fermi surface, characterized in the 
projected FL theory by a nonzero value of the traceless, 
symmetric nematic order parameter Qij = f dp Qij (p) 
where Qij{p) = {‘2-PiPj — Sij)Shp. This effectively spin¬ 
less order parameter is identical to the one that describes 
nematic order in a standard spin-degenerate FL [31]. In 
the original unprojected theory however, this translates 
into a nonzero value of Qij = J dp Qij (p) where 

Qijip) =PiSs^p+PjSsp-Sijp-6sp, (19) 

is a quadrupolar order parameter involving both spatial 
and spin degrees of freedom that was recently discussed 
in the context of possible instabilities of surface Majorana 
fermions in the topological superfluid ^He-B [32] and 3D 
spin-orbit coupled metals [17, 33]. Thus the quadrupo¬ 
lar distortion of a helical Fermi surface is necessarily ac¬ 
companied by a time-reversal invariant form of magnetic 
order similar in spirit to spin nematic order [34]. 

Nonequilibrium properties of the helical FL such as col¬ 
lective modes can also be studied using the projected FL 
theory, assuming that the relaxation-time approximation 
is valid such that scattering between states of different 
helicities can be neglected. In the hydrodynamic regime 
LOT <C 1 where r is the quasiparticle collision time, the 
helical FL supports ordinary sound waves (first sound) 
with velocity [19] 

Cl = vf\J ^{1 + Fo){l + Fi), (20) 

while in the collisionless regime ojt ^ 1 a zero sound 
mode may exist under certain conditions [25]. If Fq > 0 


only is nonzero, the zero sound velocity is given in the 
limits of strong and weak interactions by [19] 


Co ~ VF\f ^Fo, Fo OO, (21) 

Co (1 + iFo"). ^0 ^ 0. (22) 


We conclude by discussing prospects for the experi¬ 
mental determination of the projected Landau parame¬ 
ters Fi. ARPES can determine pF which, via Luttinger’s 
theorem, yields the density n. Using Eq. (13), Fq could 
then be inferred from measurements of the heat capac¬ 
ity and electronic compressibility of the surface states. 
The latter can in principle be determined directly from 
the ARPES data or via single electron transistor mi¬ 
croscopy [35]. To determine Fi, one could perform a 
transient spin grating experiment [25] to generate a spin- 
density wave with momentum q and transverse ampli¬ 
tude Sq. Due to spin-momentum locking, this will in¬ 
duce a density wave at the same momentum with am¬ 
plitude Uq. The existence of an undamped sound mode 
at frequency lo = Cgq implies a relation between the two 
amplitudes [19], 


_ 1 Cs 

Uq 1 + FiVf' 


(23) 


where Cg is either ci or cq depending on whether one 
is in the hydrodynamic or collisionless regime. Using 
Eq. (20)-(22) one can extract Fi from a measurement of 
the amplitude ratio Sq juq and previous knowledge of Fq ■ 
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This supplemental material provides a detailed derivation of the results presented in the main text. In Sec. SI, we 
justify the form of the quasiparticle interaction term used in Eq. (1) of the main text. In Sec. SII and Sec. Sill, we 
constrain the form of the interaction term by symmetries and derive Eq. (5)-(9) of the main text. In Sec. SIV, we use 
the projected Landau functional to derive the equilibrium properties of the helical Fermi liquid, i.e., Eq. (13)-(15), 
(18), and (19) of the main text. In Sec. SV we study the spin susceptibility of the noninteracting Dirac cone and 
derive Eq. (16) of the main text. Finally, in Sec. SVI we study the collective modes (sound modes) of the helical 
Fermi liquid and derive Eq. (20)-(23) of the main text. 


SI. REDUCED TWO-BODY INTERACTION 


In this section we derive the reduced two-body interaction [Eq. (1) in the main text] from a generic, translationally 
invariant two-body interaction. In this supplemental material we use the language of second-quantized interaction 
Hamiltonians, but the same reasoning applies to the interaction term in the Landau functional. A generic transla¬ 
tionally invariant interaction is given by 


V = 



cPp' 

WF 




(27r)2 Li/3;75(Pj R i ^)^p+q,a^pS^p' — S' 


(SI) 


In Fermi liquid theory, one considers the low-temperature limit where all fermionic momenta must lie on the Fermi 
surface,^ which implies that there are only three possible interaction channels: forward scattering with q = 0, exchange 
scattering with q = p' —p, and the Bardeen-Cooper-Schrieffer (BCS) channel with p' = —p but q otherwise arbitrary. 
Since we are not interested in pairing and superconductivity in the present work, we forget about the BCS channel. 
Therefore at low energies we can simplify the interaction to include only forward and exchange scattering, 


w 1 [ f 

^ {2nyJ (2^)2 


(Vaf)-^S{P,P\ 0)c],„Cp^C^,^Cp,5 -f Vcf3-jsiP,p',p' “ P)Cp'aS/3 


^PIMP'S 


/ (^ / (^“'3;75(P,P',0)cJ,„Cp^C^,^Cp,5-I4/3;75(R,P',P'“P)47Cp^4'aV5) ’ (^2) 


(27r)^ J (27r)^ 

where we have ignored one-body terms. One can show from Eq. (SI) that Fermi statistics implies 


Va/3',^y5iP,P',q) = -V^l3',a5iP,P',P' -P- q), 


(S3) 


which upon setting q = 0 yields Vap-js{p,p',0) = —V'yj 3 ',a 5 {p-,p',p' — p), i-e., the forward scattering and exchange 
scattering contributions are related. Substituting this into Eq. (S2), we find that the two contributions are in fact 
equal and simply add, 



d^p' 

i2nr 


Val3-js{P,p',0)cJ - 


'poi ^p' 7 ^p' 5 




d?p f d^p' 


(27r)^ J (27r)" 


Vap-,js{p,p')c^ ^ 


poc ^p^ ^p' 7 ^p' 5 ’ 


(S4) 


where the factor of two has been absorbed in a reduced interaction Vctp--ys{PiP') that only depends on two momenta 
p and p'. Finally, since in Fermi liquid theory we focus on momenta near the Fermi surface, we can neglect the 
dependence of Vai 3 -,-ys{p,p') on the magnitudes |p| and |p'|. One thus sets the interaction equal to its value on the 
(circular) Fermi surface, 


V ^ 



d^p' 

(27r)2 


Va/3;js(P,P)c^ '' 


pa ^pj3 ^p' 7 ^p' 5 ’ 


(S5) 


where Vap-,-yS {p,p) = Uc/3 ■.•ysiPFpiPFp') with p,p' unit vectors in the direction oip,p', andp^:’ is the Fermi momentum. 
Since the quasiparticle matrix distribution function has the same symmetry properties as the expectation value 
(CpaCp/j), the form of the Landau functional (1) in the main text follows. 
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SII. CONSTRAINING THE INTERACTION BY SYMMETRIES 


This section explains how to work out the most generic form of Vap-^s{p,p) consistent with the symmetries of the 
problem. 

Even before considering symmetries, particle indistinguishability (i.e., fermionic or bosonic statistics) gives us 

VjS-,apip',P) = Va/3-,j5{P,p'), (S6) 


and Hermiticity of V gives us 


VU,s-yiP’P') = ^a0-,js{p,p)- (S7) 

For each p,p', Va/ 3 -^js{p,p') forms a 4 x 4 matrix in spin indices. A natural basis on which we can expand this matrix 
is given by the tensor product of two sets of Pauli matrices p = 0,1,2,3 where 1 denotes the 2x2 identity 

matrix. We can write 

Va0-,js{p,p) = (S8) 

fiy 


such that in virtue of Eq. (S5), Vqo corresponds to charge-charge interactions, Vjj to spin-spin interactions, and Vq/, E/o 
to spin-charge interactions. We will use uppercase indices I, J for all three components x, y, z of spin and lowercase 
indices i^j for the in-plane components x,y. In this new basis, particle indistinguishability requires V^y{p,p) = 
V,y^{p\p) and Hermiticity requires V^u{p,p) = i.e., the matrix is real. Since the charge density and 

spin density are even and odd under time reversal, respectively, time-reversal symmetry requires 

Vbo(-p,-p ) = Voo{p,p), Vij{-p,-p) = Vij{p,p), 

for charge-charge and spin-spin interactions, and 

Voii-p,-p) = -Voi{p,p), Vioi-p,-p) = -VioiP,p), 

for spin-charge interactions. Note that Voi{p,p') = Vio{p',p) from particle indistinguishability. 

The interaction V^^{p,p') can be expanded in angular momentum components, 

OO 

V^Ap,p)= E (Sll) 

1,1'— — OO 


(S9) 

(SIO) 


where p = {cos 6p, sin 9p). Particle indistinguishability implies = Ej^* and Hermiticity implies V^J’ * = 
For charge-charge and spin-spin interactions, time-reversal symmetry implies 


J/ll' _ / 1 J/ll' 

^00 — ^00 ’ 


yll’ ^ 


J/ll' 


IJ- 


(S12) 


In other words, I +1' must be even for these coefficients to be nonzero, which is the same as saying that I and V must 
have the same parity. Therefore I' = I + 2m, m G and we can write 


'x/ll' _ Trl,l-\-2m 

^00 — *^00 


j/ll' _ Trl,l-\-2m 

^IJ — ^IJ 


m G’L. 


For spin-charge interactions, time-reversal symmetry implies 


< = -(- 1 )' 


^rLi 
^ 0 / 7 


^/o = -(- 1 ) 


y/ll' 

^10 5 


(S13) 


(S14) 


thus for these coefficients to be nonzero I + I' must be odd, which is equivalent to saying that I and I' must have 
opposite parity. Therefore V = I + 2m -I- 1, m G Z, and we have 


j/ll' _ Trl,l-\-2m-\-l _ T7'Z,/+2m+l 

^01 — *^ 0 / ’ — ^70 ’ 


(S15) 


We now turn to implementing SO{2) rotation symmetry, which is more subtle. From Eq. (S5) and (S8) we can 
write 




(S16) 
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where no{p) = n{p) is the charge density and ni{p) = 2s/(p), I = x,y, z is (twice) the spin density. Using Vqi{p, p) = 
V/o {p' I p) from particle indistinguishability, we have 


where 



V = 

Vcc + f4c + 14s ) 

(S17) 

Wc ~ 2 j 

f <Pp 

(27r)^ 

/ (27r)2^oo(p,p)^^(p)?^(p'), 

(S18) 

Vsc = 2 j 

f d'^P 

(27r)2 j 

I Mp)siip'), 

(S19) 

Vs, = 2j 

f d^p 

(27r)2 J 

f ( 2 ;r) 2 ^"*'^’^ )si{p)sj{p'), 

(S20) 


are the charge-charge, spin-charge, and spin-spin interactions, respectively. Using Eq. (Sll) and the constraints from 
time-reversal symmetry (S13) and (S15), we have 


Im ^ ^ 

Im ^ ^ 


Im 


To implement SO(2) symmetry, we note that the annihilation operator transforms as 

\ / eta' ^ ^ 


(521) 

(522) 

(523) 

(524) 


where on the left-hand side is the rotation operator, and on the right-hand side is the usual 2x2 rotation 

matrix 


R 


V — 


cos if — sin ip 
sin p cos p 


(S25) 


This implies that the charge and spin densities transform as 

R{p)n{p)R{p)~^ = n{R~^p), 

R{p)siip)R{p)-^ = R^^Sj{R-^p), for i,j = x,y, 

R{p)szip)Rip)~^ = SziR~^p). (S26) 

We require that R{p)VR{p)~^ = V. From Eq. (S26) it is clear that the charge-charge, spin-charge, and spin-spin 
interaction terms in Eq. (S21), (S22), and (S23) will transform into themselves under SO{2) rotations, and we can 
look at each term in turn. 


SI. Charge-charge interaction 

For the charge-charge interaction, requiring R{p)VccR{p)~^ = Wc implies 


E' 

Im 




7 ( 2^)2 7 ( 2 ^) 2 ® '^00 ri(pMP) 



^ P plST,+i(l+2m)9„f 
(27r)2 


T xZ ^l-\-2m 


n 


(p)n(p'). 


(S27) 
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for arbitrary (p, which implies the constraint I = —m. We therefore obtain 


Wc 



d^p' 

{2ny 


Amd. 


vv' 


00 


"n(p)n(p'), 


(S28) 


where 9ppi = 6pi — Op is the angle between p and p'. The charge-charge interaction matrix element Voo(p,p^) can 
only depend on this relative angle, i.e., 

Voo{p,p) = VooiOppi) = Y, (S29) 


In fact, because of particle indistinguishability this matrix element is symmetric under p p' and therefore depends 
only on the cosine of Oppi , 


VooiOppi) = I/p™ + 2 ^ cosmOppi. 

m—1 

Therefore the charge-charge interaction is 

= ^ j J + 2 £ V^oo™’™ cosm6lpp,^ n(p)n(p'), 


(S30) 


(S31) 


where the I/qq™’"*, 


m = 0,1, 2 ,... are real coefficients. 


S2. Spin-charge interaction 


For the spin-charge interaction, it is convenient to first separate the x, y components of spin from the 2 component, 
as they have different transformation properties under rotations, 


f4c = 2 ^ 


d^p f d?p' 




(27r)^ J (27r)^ 


{vi- 




(S32) 


Requiring Riip)VscRi‘p) ^ = I4c implies 


E' 

Im 


il2{l+i 


n.) + l]tp f ^ P f ^ P ^iW-o+i(l+2m+l)e^i 


(27r)2 J (27r)2 




'Rilnip)sj{p') + 


= E 


d^p f d^p' 


Im 


(27r)2 J (27r)2 




>' (vif+^^+\ip)s,ip') + Vt'+^"^+\ip)s.ip')) . (S33) 


For the Oz component of the interaction, this implies _ g 2 [ 2 (i-i-TO)-i-i]v 3 p^M-i- 2 m+i arbitrary p, which would 

require 2 (Z-|-to)-|- 1 = 0. This is impossible since I, m are integers, hence — q Pqj- Qa; and Oy components, 

we have 


2 [ 2 (Z+m) + l] 9 ?Ty'^,^+ 2 m+l _ -Trl,l-\-2in-\-l 


which is equivalent to 


-Qij-\A.,l-\-2m-\-\ 2 [ 2 (Z+m) + l]i;/3T7'^,Z+2m+l 


(S34) 


(S35) 


In other words, be an eigenvector of the rotation matrix R^p with eigenvalue e 42 (i-i-m)-i-i](^_ 

eigenvalues and eigenvectors of the rotation matrix (S25) are given by 




1 ( ±i 

72U 


(S36) 
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respectively. For the eigenvalue this implies that 2(1 + m) + 1 = 1 which is satisfied by I = —m. The form of the 
corresponding eigenvector implies that 


TT-—m,m+l -T/—m,m+l 


(S37) 


For the eigenvalue e this implies that 2(l + m) + l = —1 which is satisfied by I = —m — 1. The form of the 
corresponding eigenvector implies that 

T7- —m—l.m ■-tr—m—l.m /ooo\ 

K)y =*^0x ■ (S38) 

Considering these two possible values of I, the spin-charge interaction becomes 


Kc = 2 ^ 


d'^p f d^p' 


e*™®pp' 




Ore 


(27r)2 J (27r)2 

= ,.ip') + H.C. (S39) 

Unlike in a standard Fermi liquid, here there exists a time-reversal and rotationally invariant spin-charge interaction. 


S3. Spin-spin interaction 


Defining a 3D rotation matrix 


7 ? _ [ Rip 0 

/Vcip — 1 Q ^ 


(S40) 


and denoting the 3x3 spin-spin interaction matrix by y*l-i- 2 m, 3-component spin vector si(p) by 

s(p), we can write the spin-spin interaction term (S23) as 




By virtue of Eq. (S26), s{p) transforms under rotations as 


R(ip)s(p)R(ip) = n^s(R^ p). 


(S41) 


(S42) 


Therefore, requiring R(ip)VssR{p) ^ = 14s implies 

d'^p f d?p' 


E' 

Im 




■p-\-i(J,-\-2TTL)G. 


= E 


(27r)2 J (27r)2 

d^p f d^p' 


(27r)2 J (27r)2 


il9j,+i{l+2m)e^, g^p-^TYl,l+2m 


This implies = y'.'+2”^^ or equivalently 

for arbitrary p. Separating the 3x3 matrix into in-plane (x, y) and 2 components. 


Yl,l+2.m _ j || 




Tp;-/,/+2m I ’ 


(S43) 


(S44) 


(S45) 


where yj|’^+2™ is a 2 x 2 matrix, yj|>D2"‘ is a 2 x 1 column vector, and yD+ 2 m jg 7 x 2 row vector, condition (S44) 
translates into the four conditions 


p T/^i^+2mp—1 _ 22(Z+m)(/?T/'^i^+2m 

y II t V II , 

p irZ,Z+2m _ 22(Z+m)c/?T7’Z,Z+2m 

\\,z ^ \\,z ’ 

1 rZ,Z+2m o—1 _ 22 (Z+m)(/?T/•Z,Z+2m 

^ z,\\ ^ z,\\ 

^l,l-\-2m _ Z2(Z+m)c/?-^Z,Z+2m 

' r r ^ y y V 


(546) 

(547) 

(548) 

(549) 
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Beginning with the simplest condition, Eq. (S49) requires that I = —to, hence the zz part of the spin-spin interaction 
becomes 



d^p' 

{2ny 


Y-m,m 


Sz{p)Sz{p'). 


(S50) 


Similar to Eq. (S35), Eq. (S47) requires that i^g eigenvector of the 2D rotation matrix with eigenvalue 

^ As seen before, the eigenvalues of R^p are e^'’^. Since /, to are integers, condition (S47) can only be satisfied 
if = 0, thus = 0. Likewise, Eq. (S48) is equivalent to the eigenvalue condition 

(S51) 


which can only be satisfied if = 0 , thus = 0 . 

Finally, Eq. (S46) can also be converted to an eigenvalue condition by expanding the 2x2 matrix qj^ 

basis of Pauli matrices plus the identity matrix. 


irZ,Z+2m _ -%rl.l-\-27n ol 

^ II ~ ^ 


(S52) 


with a = 0,1,2, 3. The rotation matrix R^ can be written as R^, = a^cosp — icr^ sin <p, which gives the following 
transformation properties for the Pauli matrices, 

R^a°‘R~^ = cr“ (cos^ tp + (—1)“ sin^ tp) — 2(1 — 5ao)e°‘‘^^cr^ sin cos p. (S53) 

Using this property, Eq. (S46) is equivalent to the following condition. 


Tx/,Z+2m 0 I 

Eq cr + 

_ ^i2(l+m)ip 


Vi cos 2 (/?-I-Eg’ sm 


in2(/5^ 


1 I T7'/,Z+2m 2 

(7 -j- Vn (7 


(V 3 ’ 


COS 2(p— Vi sin 


in 2 (/j^ 


(^" 0 ’ 


1-^-2771 0 I T 7"^,/“|“2?71 1 I T rl ,l-\-2Tn 2 I T 7"^,/“|“2?71 3 

G V -y G -T V 2 G -|- G 


The (T° and terms give the conditions 


-rrl,l-\-2m _ tV— m,m r 

Kq — ^0 


7+m,0 5 


-ryl,l-\-2m -Tr—m^mc 
^2 — ^2 ^Z+m,05 


while the and terms give the condition 


R 


2 c/? 


l^l-\-2m 

I 


= e 


i2{l-\-m)ip 


l,l+2m 

3 

l,l-\-2m 


(554) 

(555) 


(S56) 


This is again an eigenvalue condition, but this time for the rotation matrix with angle 2p which has eigenvalues 
This implies that l + m = ±1, or I = —to± 1. The eigenvectors, however, are the same as in Eq. (S36). For I = —m + 1, 
we have 


V ~'— m+1, m+1 -tV— m+1, m+1 

3 — iVi , 

while for I = —to — 1 , we have 

T+—m—l,m—1 
^3 

The complete spin-spin interaction is therefore given by 
d^p f d?p' 


= -+■ 


—m—l,m—1 


(557) 

(558) 




which, upon substituting Eq. (S55), (S57), and (S58), yields 


= 2 ^ 


d^P f d^p j^g 


(27r)2 J (27r)2 


e pp 


+ Sys'y) + - ^.4) 

+ (e*»ps_) (e®Vs'_) - (e-®«ps+) (e-®Vs4)], (S60) 


where we denote s/ = s/(p) and s'j = si{p'). The first two terms correspond to an XXZ interaction, the third term to 
a Dzyaloshinskii-Moriya interaction which reflects the presence of spin-orbit coupling in the system, and the last two 
terms to anisotropic spin-spin interactions similar to those found in compass models,^ but with a continuous rather 
than discrete symmetry. These terms are not invariant under separate spatial and spin rotations, but only under a 
simultaneous rotation in spin space and real space. 
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Sill. LANDAU PARAMETERS 


Landau parameters are the real coefficients of Hermitian interaction terms. The charge-charge interaction (S31) is 
already in this form, but we wish to write the spin-charge interaction (S39) and the spin-spin interaction (S60) in this 
form as well. We first consider the spin-charge interaction. Splitting into real and imaginary parts, 


T 7-—?7i,m+l 

^Ox 


= v' + iV" 

^ m. ' ^ m 


(S61) 


Vac can be written as 




d^p f cPp' 


(27r)2 J (27r)2 


-h 


OO 

^0 + cosmOpp, - (U" - U"„) sinm 6 »pp/) 

m—1 

oo 

Vq -t- ^ [(U" -H U"„) cosmOpp, + {V^ - v!_^) sinmOpp 


n{p)p ■ s' 


n{p)pxs'), (S62) 


where we have used the fact that = Px ± ipy, and we denote z ■ {a x b) = a x b for simplicity. For the spin-spin 

interaction, particle indistinguishability implies that Uq”’"’™, 14 “™’™, and are real, while = 

Splitting ]-gai ajjfj imaginary parts. 


^_m+l,m+l ^ ^yrny 


(S63) 


the spin-spin interaction can be written as 


Ks = 2 


d^p f cPp' 


(27r)2 J (27r)2 


Uq™ + 2 ^ I 4 cosmOppi j s • s' + ( 14 °° + 2 ^ cosmOppi j s^s'^ 


m—1 


m—1 




2 " sin mOpp! s x s' 


+ 2 1 ^( 14 °)' + 2 ^ (V]™)' cosm6»ppj ((p • s)(p x s') + {p x s){p' ■ s')) 
- 2 ( ( 14 °)" + 2 ^ ( 14 ™)" cosm 6 »pp/ ] ((p • s)(p • s') - (p x s)(p x s')) 


(S64) 


Considering the full interaction term V = 14c + 14c + 14s; for each m there are ten independent real coefficients, hence 
ten Landau parameters. We define one charge-charge Landau parameter /“; 


TCC _ 

dm 

four spin-charge Landau parameters ..., 

xsc,l _ 

J m 

rsc,2 _ 

J m 

tsc,3 _ 

J m 

tsc,4 _ 

J m 




m = 0 


2144 ’"’™, m = 1,2,3,... 


414 , m = 0 

+ m = 1,2,3,... 

0 , m = 0 

-4{VZ-V!!^), m = 1,2,3,... 


(S65) 


414", 


m = 0 


4(14( + 14'„), m = 1 , 2 , 3 ,... 

0, m = 0 

4(K-Vlm), m = 1,2,3,... 


(566) 

(567) 

(568) 

(569) 












and five spin-spin Landau parameters ■ ■ ■, 


?ss,l _ 

r dEgOO, 

m = 0 

m 

\ 

m = I, 2, 3,. 

7SS,2 _ 

1 4H™, 

0 

II 

m 


m = 1, 2, 3,. 

7SS,3 _ 

/ .0, 

TO = 0 

m 

\ 8 ^ 2 ”'"’™, 

TO = 1, 2,3,. 

?ss,4 _ 

j 8 ( 1 / 0 )', 

TO = 0 

m 

1 i6(vry, 

TO = 1, 2, 3,. 

?ss,5 _ 

f -8(ViT 

, TO = 0 


[ -W{Vry\ m = l,2,3,... 

In terms of these Landau parameters, the interaction terms can be written as 


= i E 


d^p f cPp' 


m—0 '' 


(27r)2 J (27r)2 


/“ cos m0ppm{p)n{p') 
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(570) 

(571) 

(572) 

(573) 

(574) 

(575) 


for the charge-charge interaction, 

OO 


m—0 


f <Pp 

f <Pp' r 

j (27r)2j 

(2^)2 L 


(/m’^ cosmOpp, + sinm6»pp/) n{p)p ■ s{p') 

+ (/m’^cosm6»pp/ -h sinmepp>) n{p)p x s[p') , 


(S76) 


for the spin-charge interaction, and 

- r ^ r dy 
(27r)2y (27r)2 

+ /m cosmOpp, [(p • s{p)) {p X s{p')) +{px s{p)) (p ■ s(p'))] 

+ /m cosmOpp, [{p ■ s{p)) {p ■ s{p')) - (p X s(p)) (p X s(p'))] }, (S77) 

for the spin-spin interaction. Interpreted as quasiparticle interaction terms in a Landau functional, Eq. (S75)-(S77) 
correspond to Eq. (5)-(7) of the main text. 


{/m cosm0pp:s{p) ■ s{p') + cosm0pp's^{p)s^{p’) + smm0pp,s{p) x s(p') 


SI. Projected Landau parameters 

In this section we explain how to derive the projected Landau functional [Eq. (8) of the main text] from the 
unprojected theory we have just described. The starting point is to drop the negative helicity part in the expressions 
for the fermion operators, 


le 


Cp-f — 


V2 


" (V'p+ V’p-) 


le 


iO-n 


V2 


■'0P+) 


Cpi — V'p-) ~ 


(578) 

(579) 


which can be expressed as Cp„ ~ ’rjp^tjjp where the c-number spinor pp = -^(*6 1) obeys PpPp = I, and we define 

the effectively spinless fermion operator ^j;p = ipp+. One then substitutes this expression for Cpa into the interaction 
Hamiltonian (S5). This produces a projected interaction Hamiltonian V, 




/ 





(27ry 


^(P,P)V’p'0pV'p'^p', 


(S80) 
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where the projected matrix element V{p,p) is given in terms of the unprojected ones by 

V{p,p) = Va^.,psip,p)vpaVpf3Vp'jVp's = P)iVp<^^Vp){Vp'. 


(S81) 


where we used Eq. (S8). The quantity Vp^^Vp can be thought of as the expectation value of in the single-particle 
eigenstate at p on the Fermi surface. We have ? 7 |,cr°r;^ = 1 which corresponds to a particle number of one, Vp^^Vp = 0 

which indicates that spin polarization on the Fermi surface is entirely in-plane, and Tip<J^Tip = tijPj, i,j = 1,2, which 
indicates that spin is perpendicular to momentum everywhere on the Fermi surface, i.e., spin-momentum locking. 
Using the fact that Voi{p,p) = Uo(p^p), and the fact that the interaction (S80) must be symmetric under p o p, 
we have 


V(p, p) = Vbo(P, P) + 2Vo^(p, P (P> p)(^ikejiPkPi■ (S82) 

The Vzz component of the spin-spin interaction does not enter since spins on the Fermi surface are entirely in-plane. 
We can now read off Voo{p,p), Vbi(p,p^), and Vij{p,p) from the interaction terms (S75), (S76), and (S77), and their 
original definitions (S18), (S19), and (S20), respectively. We find 


vooiP, 


m—0 
oo 


^ oo 

Voi{p,p ) = 2 ^ cosmOpp, -f sinmOpp') p\ 

m—0 
oo 

yj{p,p) = 4 X! cosmepp,Sij + sinm0pp>e^J 

m—0 


m—0 

T /m’ cosmQpp! {ekjPiPk 

Eq. (S82), we 


“C Jm vun//«/pp' 

Substituting these expressions into 


CkiPkPj) 


'^'j) + /m’^ cosmOpp, {piPj - ekiCijPkPi) ■ 

find 


(583) 

(584) 


(S85) 


^(P, P ) = X] f’- ^^PP' ’ 

1=0 

where the projected Landau parameters fi are given by Eq. (9) in the main text. 


SIV. EQUILIBRIUM PROPERTIES OF HELICAL FERMI LIQUIDS 


This section presents a detailed derivation of the physical properties of helical Fermi liquids from the projected 
Landau functional [Eq. (8) of the main text]. Rather than as the coefficients of terms in a second-quantized interaction 
Hamiltonian operator, we would really like to think of fi as the coefficients in this functional, 

-^)(5hp + -Xy J fi cos Wpp^Snpdfip,, (S87) 

where = vpp and we have explicitly added a chemical potential term p. The (renormalized) quasiparticle energy ep 
is given by the functional derivative of the Landau functional with respect to the quasiparticle distribution function. 


_ 6E 

^ drip 



cPp 

(2^)2 


fi cos lOppf 


Sn. 


p' • 


(S88) 


It is important to note that vp is the renormalized Fermi velocity and not the bare (noninteracting) one, which is 
denoted by v^p. 
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SI. Specific heat 


We first investigate the specific heat. The derivation we use closely follows the standard derivation for the specific 
heat in standard Fermi liquids.^ The entropy density, s, is given by 


J (2^ - (1 - «p) ln(l - rip)), 


s = -ks 


where ks is Boltzmann’s constant. The variation is the entropy density is given by 


Ss = —ks 


-^dn In 
(27r)2'’ P 


1 — Ur 


1 


(Pp 

tJ ^ 


6np{ep - p), 


where T is the temperature. The variation in particle density can be written as 

Cp- 6^ , 


c— Oflp 


T 


-ST + Sep — p . 


To lowest order in T, we have 


Ss = - 


1 f cPp dup 


(cp - p)^ST = -klg{p) 


de 


e — dn TT 


T2 J (27r)2 dcp 

where p(p) is the density of states at the Fermi surface which is given by 

d?p 


ksT J 


de 




P(.) = / 


(27r)' 


:S{e-ep). 


For the noninteracting case we have p{p) = p,/2Trvp. The specific heat is then 

TT^ 

Cv = -jp{p)k%T. 

One then defines the electronic specific heat coefficient 7 as the zero-temperature limit of CyjT^ 

7 = \rr‘^kpBP{^F)- 


(S89) 


(S90) 


(S91) 


(S92) 


(S93) 


(S94) 


(S95) 


S2. Compressibility 


We now turn to the electronic compressibility. Again, this derivation closely follows that for standard Fermi liquids. 
The compressibility k at zero temperature is defined as 


1 dn 

rp dp' 


(S96) 


where n is the density of electrons, given hy n = f Pp where Pp = ^^^(‘"jxTCpcr)- Projecting the field operators 
to the Fermi surface, we obtain pp = rip as expected, thus Spp = Sfip. At zero temperature, the density variation 
[Eq. (S91)] is 


Srip = ^^{Sep - Sp). 

0€p 


(S97) 


The quantity vanishes everywhere except at the Fermi surface and the variation of p produces a variation of Srip 


de 


that is isotropic. Integrating Eq. (S97) over momentum, we find 

dPp dfi^ 


Sn = 


X- _ 


d'^p drip 
(27r)2 dep 


{Sep - Sp) = 


(27r)2 dep 


E 

' 


d^p' 
{2Tr)' 


~fi cos lOpp/Sripf Sp J . 


(S98) 
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For the last step, we used Eq. (S88) for Scp. After integration over p', only the I = 0 contribution remains and we 
have 


cPp dup I f cPp' - _ ,n s- r 1 f 


E 


J (27r)2 dep \f^J (27r) 

Defining dimensionless projected Landau parameters Fi as 


fI cos I0pp> Slip/ Sfjj j — 


(27r)2 de^ 


{fo6n - 5p). 


Fi = p{eF)fi 


de 

27T 


COS W = < I 


p(ep’)/o, 1 = 0, 


y{eF)fi, Z = 1,2,3,... ’ 


we arrive at 


) 1 _ 

n2 1 + Fn 


(S99) 


(SlOO) 


(SlOl) 


S3. Spin susceptibility 


We now investigate the spin susceptibility of a helical Fermi liquid. As mentioned in the main text, contrary to 
a standard Fermi liquid here the spin susceptibility is not strictly speaking a Fermi surface property. Therefore the 
present projected Fermi liquid theory can only correctly describe the spin susceptibility of the helical Fermi liquid in 
a certain limit to be seen below. ^ 

The total spin density of the helical Fermi liquid is s = f Sp where 

«P = ^(CpO-Cp)- (S102) 

Projecting the fermion operators to the Fermi surface, we obtain 

■Sp = ^e*jTjAp, * = 1,2, (S103) 

s^p = 0. (S104) 

Therefore, we have 6Sp = ^eijPjSfip and Ssp = 0. As a result, our projected Fermi liquid theory will predict a zero 

out-of-plane susceptibility Xzz = 0. For the in-plane susceptibility, consider applying an in-plane Zeeman term, 

6H = - J ■ ic^crcp, (S105) 

where pb is the Bohr magneton and g is the g-factor of the helical Fermi liquid. Ignoring constant terms, this leads 
to a change in the energy. 


6E = {6H) = - 

where the change in quasiparticle energy is 


d?p 


ggsB ■ 5sp = 


d?‘p 

IF? 


Scpi^B^Sfip^ 


Sep(B) = ep(B) - Cp = --gfinB x p. 


(S106) 


(S107) 


The variation of Sup is given by 


Fin 

Sfip = -^{6ep - 6p), 


(S108) 


where 


1 


OO 

<5ep = -^gPsBieijPj + ^ 


d^p' j 


(2^)^ 


fi COS lOppfSfip/. 


(S109) 
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Since the chemical potential p is a scalar and does not depend on the direction of the magnetic field, its variation can 
be ignored when calculating the linear susceptibility. We now introduce a renormalized g-factor gi{p) that depends 
on quasiparticle momentum 


1 1 oo 

Sep = = --gpsB^eijPj + ^ 


d^p' I 


1=0 • 


{2ny‘ 


fi cos lOppf Svipf . 


Inserting Eq. (S108) into Eq. (SI 10), we find an integral equation for gi{p), 


9i{p) 


9 


= 


d^p' j dup, g,{p') 

; Jl cos Wpp'- 


(SllO) 


(Sill) 


i=0 ' 


(2^)2-— g 

We note that integral equations also appear for a partially spin-polarized Eermi liquid.^ The spin susceptibility is then 
9PB 


Xi^ = lim 


B->0 B 


'i'!’ = lA 


(27r 


(Pp ^ dup I ^ 


(2^) 


1=0 ' 


d^p' j dup, 

fI cos I9pp,^^g^{p) 


dcp! 


To make some progess in the interacting case, we assume a solution of the form 

pipy — PeHeikPk ■ 


(S112) 


(S113) 


Only the I = 1 term will survive. The I = 1 cosine term can be written as cosOpp' = PxP'x + PyP'y After evaluating 
the angular integral in Eq. (Sill), we find 


Heff — 


9 


l + Fi 

Turning to the spin susceptibility given by Eq. (S112), we find 


5 VI 


" 4(1 + Fi) y (2^)2 dcp " 8^ PBPiep) ^ ^ , 


(S114) 


(S115) 


at zero temperature. We also find that Xij vanishes for i ^ j, which can be explicity seen from Eq. (S112). This can 
be compared to the spin susceptibility of the noninteracting helical Eermi gas, which is derived using both helicities 
in Sec. SV. 


S4. Pomeranchuk instabilities 

In this section we investigate the stability of the Fermi surface. The distortion of the Fermi surface can be charac¬ 
terized by an angular dependent Fermi wavevector,^ 


00 

PF{e) - PF = Y. ■ (S116) 

I—— 00 

The change in energy is then 

00 

SE[5np]="^Y^^ + P^)\M\ 

^ 1=0 


(S117) 


Here we have used that fact A* = A_/ since pf{9) is real. The Fermi surface is stable against spontaneous distortions 
only if 5E > 0, i.e., if Fi > —1 for all 1. Let us discuss briefly the special case of the I = 2 Pomeranchuk instability, 
which corresponds to a nematic instability.® The I = 2 projected interaction is of the form 



dPp 


(fp' 

(27r)2 


cos 20 ppf STipSiT/p^ — 


4 


(pp 


-^^TrQ{p)Q{p'), 


(S118) 


where 


QijiP^ — y^PiPj Sij^Sflp^ 


(S119) 
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from which we can construct a traceless, symmetric 2D nematic order parameter® Qij = J Interestingly, 

Eq. (9) in the main text shows that this type of interaction can be obtained from p-wave unprojected spin-spin 
interactions, i.e., the unprojected Landau parameters and . In fact, as mentioned in the main text [see 
Eq. (19)] one can construct a 2D nematic order parameter in terms of the unprojected spin degrees of freedom, 

Qij{p) = PiSp + PjSsp - S^jp ■ dsp. (S 120 ) 

This type of nematic order parameter was first considered in Ref. 6 as a possible instability of 2D Majorana fermions, 
and its 3D analog was considered in Ref. 7 in the context of spin-orbit coupled 3D metals. One can then show that 


1 f d?p 

2] 


(fp' 

(2^)2 


TrQ{p)Q{p') = 


cPp f cPp' 


(27r)2 J (27r)2 


(^cosOppfSsp ' Ssp' sin0pp'Ssp X Ssp/^. 


(S121) 


which, comparing with Eq. (S77), corresponds to a spin-spin interaction with = —/i*’^ 7 ^ 0. Because of Eq. (9) 
in the main text, this corresponds indeed to a nonzero contribution to / 2 . In fact, if we project Qij{p) to the Fermi 
surface in the sense of replacing dSp by its expectation value on the Fermi surface (Ssp) = ^CijPj, we obtain 


iQ^jip)) 


1 / 2p^Py pI-PI\ 

2 \fy-pl -‘2pxPy J ’ 


(SI22) 


which is essentially equivalent to Eq. (SI 19) except for a rotation by 7 r /4 about the z axis: by rotating p —>■ Rt^/^p, 
we have {Qij{p)) \[‘^PiPo - 5ij). 


S5. Renormalized velocity 


In this section we consider the renormalization of the Fermi velocity. We begin by considering a microscopic 
Hamiltonian in first quantization, 

H = VpZ ■ {a- X p) + Hint, (SI23) 

where the electron-electron interaction H^t is assumed to not depend on momentum. The renormalization of the 
Fermi velocity is similar, in spirit, to the renormalization of the quasiparticle mass in normal Fermi liquids. The 
derivation, however, is quite different because spin-orbit coupling breaks Galilean invariance. Following Ref. 8 , we use 
the fact that the total flux of quasiparticles is equal to the total flux of particles. To hud the velocity operator of the 
particles, we use the commutation relation 

Vf, = —^[x,H] = Vp{z X (t). (SI24) 

Because the interaction is momentum-independent, we have [x^Hint] = 0 and the velocity operator is the same as in 
the absence of interactions. By equating the total flux of particles and quasiparticles we find 

/ (^);f2 X (V'pCrV'p)) = J -^^^pVpep. (SI25) 

After projecting the fermion operators to the Fermi surface and varying both sides of Eq. (S125), we find 


/ /-/2 r ^2 p fj‘2, p ^“2. / 

j^v^ppSnp = J ^^(Vpe°)< 5 hp -J2j J cos/ 0 pp'(Vph° )^hp. (SI 26 ) 


After relabeling p p' and equating the integrands, since the variation of Up is arbitrary, we find 


VpP 


r ^ 2 p/ _ / r ^ 2 p/ _ \ 

= (2^'^' cos Wpp,Wp>n°p, = yF + ( 2 ^-^' cos Wpp^S{vFP - p)p j ■ 


Multiplying both sides by p, we obtain 


(SI27) 


/ ^ y. ^2pi _ \ _ 

= "CF 1 + ^ cosl9pp'5{vpp - p) cosOpp, ) = uf(1 + Fi). 


1=0 • 


{2nf 


(SI28) 


Only the I = 1 term contributes, and we arrive at Eq. (14) of the main text. 
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SV. SPIN SUSCEPTIBILITY OF THE NONINTERACTING HELICAL FERMI GAS 

In this section we calculate the spin susceptibility of the noninteracting helical Fermi gas while taking both helicities 
into account (i.e., without projecting out the negative helicity part). 


SI. Out-of-plane spin susceptibility 


We first consider the out-of-plane susceptibility at zero temperature. We consider a free Dirac system with a Zeeman 
term, ignoring orbital effects of the magnetic field. (The combined effects of Zeeman and orbital couplings on the spin 
susceptibility were studied in Ref. 9.) Landau quantization is expected to dominate only at very low fields. Spefically 
the orbital contribution will dominate if 


hvp 



> gfisB. 


(S129) 


Experimental parameters for 612803 , Sb 2 Te 3 , and 612 X 63 exposed to ambient conditions put this scale on the order 
of 10“^ T, thus for those systems, the Zeeman effect will dominate under typical experimental conditions.(When 
exposed to ambient conditions the Fermi velocity can decrease by two orders of magnitude, which allows the Zeeman 
term to dominate down to very small fields.) The Hamiltonian of the system is given by 


H = 


(Pk I 


{h{k) - p 


igpsBcr^) Cfc, 


where Ck = {ck^,Cki) is a two-component Dirac spinor, and 


h{k) = VpZ ■ {(T X k) = Vp 


0 ike \ 

0 ) ■ 


(8130) 


(8131) 


We note that the out-of-plane Zeeman term is not diagonal in the helicity basis and thus cannot be captured by our 
theory. This Hamiltonian can be diagonalized exactly. The full Hamiltonian can be written as 


H = v°p 


(27r)2 


Ck- 


(8132) 


The eigenenergies E^{k) = xVpk+0{B^) do not change to linear order in out-of-plane field strength. The eigenvectors 
for a given chirality y = ±1 are 


li’xik)) 


^ 1«(1 - *«ff) j ^ 


(8133) 


We now calculate the expectation value of the Pauli matrices for a given chirality to lowest order in field strength. 
This expectation value is proportional to the magnetization. We find 


= X (^1 + sindfc, 

8 umming over momentum and chirality, we find 

{(Jz) = -J2j -^;;^('^xik)Wz\^xik))nFiE^ik)) = ’ 


(8134) 


(8135) 
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where npis) = + 1) ^ is the Fermi function (evaluated at zero temperature j3 —>■ oo), and {ux) = {<Jy) = 0 

due to the angular integral vanishing. We have also introduced a high-energy cutoff A. Evaluating the integrals, we 


find 




47r(-i;^)2 


(A - ep) 


2 


[p(A) -p(eF)], 


(S136) 


where p{e) = |e|/27r(u^)^ is the density of states of the helical Fermi gas. This corresponds to an out-of-plane 
susceptibility 


Xzz = = ^g^pl[p{A) - p{eF)]. (S137) 

In the limit that ep —>■ A, the out-of-plane spin susceptibility vanishes in agreement with the projected helical Fermi 
theory. 


S2. In-plane spin susceptibility 


In this section we consider the in-plane susceptibility. We take the in-plane magnetic field to be in the x-direction 
without loss of generality due to SO{2) rotation symmetry. The Hamiltonian we consider is 


H 


/ (fk 

j^cl {h{k) - p- ypsBa^) c^. 


The energy of an eigenstate of given chirality to linear order in field strength is 

Bx{k) = xv°pk - +0{B^), 

and the eigenstate of a given chirality is 




(S138) 


(S139) 


(S140) 


We now calculate the expectation values of the Pauli matrices. We find to linear order in field strength 

{^^{k)\a^\'ipx{k)) = 0, 

= x(^sin6lfc -cos^ )’ 


{il^^{k)\ay\ijj^{k)) = x( cosOk - cosOkSinOk 


ggBB \ 
2vpk ) 


Summing over momentum and chirality we find 

cPk 

j ( smdfc - cosdfc^ 


(cTx) = J (^sindfc - 0 ^e_F - Vpk + sin6>fc ^^^^ ^ - 0 ^e_F + v°pk - 


(S141) 


(S142) 


Simplifying the equation by using the fact that 0(eF + Vpk — sin 0k ^^2^ ) is always one for a large Fermi energy ep, 


we find 


Performing the integration, we find 


which gives for the susceptibility 


rkp+sinffh 


dk 
2 p 


"W dk 


— I 1 sindfcfc —cos dfc 


ggBB \ 

2v% )■ 


gpsB 

i— py^)^ 


Xxx = -g‘^p%p{A). 


(S143) 


(S144) 


(SMS) 


This agrees with the result obtained from our projected helical Fermi liquid theory [Eq. (S115) with Fi = 0] when 
ep —>■ A. 
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SVI. COLLECTIVE MODES IN A HELICAL FERMI LIQUID 

In this section we investigate the collective modes of a helical Fermi liquid in the presence of a monochromatic 
external scalar potential, U{r,t) = with q pp and w <C ep- = vppp. The quasiparticle distribution 

function Up obeys the kinetic equation 

dfi (t* 

— t) ■ Vr.np(r, t) - Vr.ep(r, t) ■ \7pnp{r, t) = /[tip], (S146) 

where I[np\ is the collision integral. In the presence of an external scalar potential, np{r,t) and the quasiparticle 
energy are given by 

hp(r,t) =n°+(5hp(r,t), ep(r, t) = e° + [/(r, t) + ^ fi cos Wpp>6np>. (S147) 

In general, the collision integral involves scattering between states with different helicities, and one must keep both 
helicities. However, if the relaxation-time approximation is valid, scattering between states with different helicities 
can be neglected. 


SI. First sound 


We first consider collective modes in the hydrodynamic regime (regular sound waves), i.e., lot ^ I where r is the 
quasiparticle collision time. Our goal is to find the sound velocity in the presence of quasiparticle interactions. This 
derivation does not follow the standard derivation for regular Fermi liquids due to the lack of Galilean invariance. We 
first obtain the local momentum conservation law. To do so, we first multiply Eq. (SI46) by pi and then integrate 
over p, which gives 


where 


dg^{r,t) dT,j{r,t) 
dt drj 


cPp 


d€p{r,t) 

dri 


np{r,t) = 0, 


(8148) 




(Pp 

w? 


Ptnp{r,t), 


Tiji'i'i t) 


cPp 


dep{r,t) _ 
dp, 


(rp), 


(8149) 


where gi is the momentum density. The integral of the collision term vanishes due to conservation of quasiparticle 
momentum. We can rewrite this equation as 


^g^ir,t) dll,,{r,t) .dU{r,t) 

at + dr, dr, ~ 

where n{r,t) = / -^;^np{r,t) and H^ is the total stress tensor, given by 

Hij = Tij — Sij J ^2^)2 0- 

In general, the stress tensor is of the form 

Hij — PSij Oij , 


(8150) 


(8151) 


(8152) 


where P is the pressure and atj is the dissipative part of the stress tensor. We neglect a for the rest of this work 
as it does not have an effect on the sound velocity. We now turn to the local energy conservation law. Multiplying 
Eq. (8146) by ep(r,t) and then integrating over p gives 


Pp 

(271-) 


2 


,t) 


dnp{r,t) d 


dt 


drj 


Pp , .dep{r,t)_ 

Cp{r,t )—^- np{r,t) = 0. 


(27r) 


dp. 


We now linearize Eq. (8150) and (8153). First we define. 



dpp 

(27r)2 


Snp{r, t) = Sn{p, r, t) 



dpp 


pSnp{r,t) = SV,{p,r,t), 


(8153) 


(8154) 
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and mode expand Sn{p,r,t) and Sft{p,r,t) as 

OO OO 

Sn{p,r,t)= ^ , Sn{p,r,t) = ^ Cn{r,t)e"'^^. 


(S155) 


n——oo 


n— — oo 


The total density and energy fluctuations are then 

Sn{r,t) = 2TTBo{r,t), 6e{r,t) = 2TTVFCo{r,t). (S156) 

We note that the energy and density fluctuations are related by the equation of state. To linear order, we have 


3 f^( 71^1 

Seir^t) = —^-^6n{r,t) = n6n{r,t), 


(S157) 


where is the equation of state in the noninteracting limit. Upon linearizing Eq. (S150) we obtain 


-pdjn{r,t){l + Fq) + 2TTdt{ReCi{r,t)) = -n°d^U, 


(S158) 


and 


-^dy6n{r,t){l + Fq) — 27r9t(ImC'i(r, t)) = —nPdyU. (S159) 


Here we have used 5P = which is valid since the system is in local thermodynamic equilibrium, and 


djSP{r,t)6ij = ^p^^6n{r,t){l + Fo), 

which follows from Eq. (SlOl) for the compressibility. Linearizing Eq. (S153) gives 

dtS€{r,t) + 2TTVpda;{ReCi{r,t)){l + Fi) - 2TTVpdy{lmCi{r,t)){l + Fi) = 0. 


(S160) 


(S161) 


Taking the temporal derivative of Eq. (S163), the x-derivative of Eq. (S158) and the y-derivative of Eq. (S159), and 
substituting Eq. (S158) and (S159) into Eq. (S163), we obtain 


d^Se{r,t) - ^pvpW'^dn{r,t){l + Fi){l + Fq) = + Fi ). 

Using Eq. (S157), we find the equation of motion for the density fluctuations to be 

dMr, t) - t)(l + Fi)(l + Fq) = ^n°V2c/(l + Fi). 

Thus the velocity of first/hydrodynamic sound in the presence of quasiparticle interactions is given by 


Cl = vpJ\{\ + Fo)(l + Fi). 


(S162) 


(S163) 


(S164) 


S2. Zero sound 


We now turn to the collisionless regime ojt ^ 1. Solving Eq. (S146) by Fourier transform for U = 0, we find 


dn° 


{w - q ■ Vp)5np + q ■ Vp-^ ^ 

OCr 


Z=0 ' 


d^p' 

{2ny 


fi cos 0pp> Stip' — 0, 


(S165) 


where Vp = vpp. Following the standard approach in regular FLs, we assume a solution of Sup of the form 


dnl 

dep 


(S166) 
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and expand as 


= E 


Jm9- 

C t^r) 


(S167) 


We then have 


{uj-q- Vp)vp -qVpY^ p{^F)fi 


r2Tr 


de. 


1=0 


2tt 


— cosWppiVpi = 0 . 


(S168) 


Choosing a system of coordinates for p such that q ■ p = qcosOp, and dehning the dimensionless variable s = ui/vpq, 
we obtain 


cos 0TJ 


s — cos 9, 


■ E 


dSp' 

• COS 10 ryn'l^r)' = 0. 


/=0 


27T 


We first consider the case that Fq only is nonzero. Integrating over both Op and we obtain the equation 


(l + Foflo^s)) i/Q = 0, 


where we define the dimensionless function 


riiis) = 


dO cosOcoslO 


2tt cos 6 — s 

which is easily evaluated for the first few values / = 0,1, 2, 

no(s) = 1 - 


= ^-i{s), 


A nontrivial solution lyg ^ 0 requires 


which is easily solved to give 


\/s2 - 1’ 
fll(s) = sfiois), 
fi 2 (s) = l + (2s^-l)no(s). 


1 + AQfl(s) — 0, 


UJ 

vpq 


l + Fd 




Co_ 
Vf ’ 


(S169) 

(S170) 

(S171) 

(S172) 

(S173) 

(S174) 

(S175) 

(S176) 


where cq is the velocity of zero sound. A zero sound mode thus exists for all positive (repulsive) values of Fq. Because 
i2o{s) is real only for s > 1, for an undamped mode one must restrict oneself to s > 1. Simple expressions can be 
obtained in the limits of strong and weak interaction, 


Co-Up-y—, Fq^OO, 

Cq ^ Vf W / ’ -^0 0- 


(S177) 

(S178) 


We now consider turning on a nonzero value of Fi in addition to a positive Fq. One then obtains three coupled 
equations for vq,Di,D-i, 


(l + FoOo(s)) Fq + Fifli{s){Fi + F-i) — 0, 

2 ^oOi(s) 7 'o + (l + AiOo(s) + 7^1122) + f - i ) = 0, 

(1 + FiOo(s) - Fi02(s)) {fi - F-i) = 0. 


(S179) 

(S180) 

(S181) 
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We see that the I = 0 mode vq and the symmetric combination of the Z = ±1 modes vi + v-i are coupled by the 
first two equations, while the antisymmetric combination i?i — P-i decouples. The equation for the latter would also 
be found in a model with a pure Fi interaction. Since we are primarily interested in the effects of a nonzero Fi 
interaction on the I = 0 mode found earlier, we will focus on the first two equations (S179)-(S180). The condition of 
a nontrivial solution for z/q and + i^_i gives 


1 _ (1 + Fi)^o{s) 

“ ~l + Fi(l + 2s211o(s))- 


(S182) 


Consider negative values of Fi. One can show that Oo(s) and 1 + 2s^Oo(s) are negative for all s > 1. For Fi < 0, we 
have 


1 (l-|Fi|)|^2o(g)| 

Fo 1 +|Fi||l + 2s2no(s)r 


(S183) 


The right-hand side of this expression becomes negative for Fi < —1, implying that zero sound is destroyed for 
sufficiently attractive values of Fi. Given Eq. (11) in the main text, this can occur, for instance, due to sufficiently 
attractive microscopic I = 0 interactions in the spin channel, e.g., /q^’ sufficiently negative. The disappearance of 
zero sound due to sufficiently attractive interactions in the spin channel was also found in a microscopic study of the 
helical Fermi liquid. 


S3. Determining Fi from first/zero sound 


As discussed in Sec. SIVS5, there is an operator identity that relates the electron velocity operator Ve to the 
electron spin cr. This identity is valid in the presence of interactions, but only involves the noninteracting Fermi 
velocity Vp. Combined with the continuity equation dtUq = —iqjq where Uq is the density operator and jq = q ■ jq 
is the longitudinal current density operator, this yields the identity dtUq = —ivpqSq where Sq = z ■ {cr x q) is the 
transverse spin density operator. Passing to the frequency domain, this gives^^ 


T 

^ 

D ’ 

Vpq Uq 


(S184) 


where the right-hand side is now interpreted as a ratio of expectation values. In Ref. 12, the authors suggest generating 
a spin-density wave of momentum q and amplitude with a spin grating. In the presence of a collective mode of 
frequency to = Cgq where Cg is the sound velocity (ci or cq depending on whether one is in the hydrodynamic or 
collisionless regime), this will generate a long-lived density wave at momentum q whose amplitude Uq can in principle 
be measured. Using Eq. (S184), the ratio of amplitudes of the original 
should be given by 

Sq Cg 1 Cg 
nq 1 + Fi VF ' 


spin-density wave and induced density wave 

(S185) 


Assuming for instance that one is in the hydrodynamic regime lot <C 1, one would get 


' I + Fq 
^ 1 -I- Fi 


(S186) 


such that the value of Fi can be extracted from a measurement of the amplitude ratio, assuming that Fq is known 
from heat capacity and electronic compressibility measurements, as explained in the main text. While the sound 
modes give a g-independent ratio of amplitudes, the spin plasmon mode^2 ^j^g unscreened Coulomb interaction 

gives a ratio s^juq oc l/y^, which can in principle be used to discriminate between the two types of collective modes. 
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